3.1 \(\int \frac{\sin ^4(x)}{a+b \sin (x)+c \sin ^2(x)} \, dx\)

Optimal. Leaf size=323 \[ -\frac{\sqrt{2} \left (-\frac{2 a^2 c^2-4 a b^2 c+b^4}{\sqrt{b^2-4 a c}}-2 a b c+b^3\right ) \tan ^{-1}\left (\frac{\tan \left (\frac{x}{2}\right ) \left (b-\sqrt{b^2-4 a c}\right )+2 c}{\sqrt{2} \sqrt{-b \sqrt{b^2-4 a c}-2 c (a+c)+b^2}}\right )}{c^3 \sqrt{-b \sqrt{b^2-4 a c}-2 c (a+c)+b^2}}-\frac{\sqrt{2} \left (\frac{2 a^2 c^2-4 a b^2 c+b^4}{\sqrt{b^2-4 a c}}-2 a b c+b^3\right ) \tan ^{-1}\left (\frac{\tan \left (\frac{x}{2}\right ) \left (\sqrt{b^2-4 a c}+b\right )+2 c}{\sqrt{2} \sqrt{b \sqrt{b^2-4 a c}-2 c (a+c)+b^2}}\right )}{c^3 \sqrt{b \sqrt{b^2-4 a c}-2 c (a+c)+b^2}}+\frac{x \left (b^2-a c\right )}{c^3}+\frac{b \cos (x)}{c^2}+\frac{x}{2 c}-\frac{\sin (x) \cos (x)}{2 c} \]

[Out]

x/(2*c) + ((b^2 - a*c)*x)/c^3 - (Sqrt[2]*(b^3 - 2*a*b*c - (b^4 - 4*a*b^2*c + 2*a^2*c^2)/Sqrt[b^2 - 4*a*c])*Arc
Tan[(2*c + (b - Sqrt[b^2 - 4*a*c])*Tan[x/2])/(Sqrt[2]*Sqrt[b^2 - 2*c*(a + c) - b*Sqrt[b^2 - 4*a*c]])])/(c^3*Sq
rt[b^2 - 2*c*(a + c) - b*Sqrt[b^2 - 4*a*c]]) - (Sqrt[2]*(b^3 - 2*a*b*c + (b^4 - 4*a*b^2*c + 2*a^2*c^2)/Sqrt[b^
2 - 4*a*c])*ArcTan[(2*c + (b + Sqrt[b^2 - 4*a*c])*Tan[x/2])/(Sqrt[2]*Sqrt[b^2 - 2*c*(a + c) + b*Sqrt[b^2 - 4*a
*c]])])/(c^3*Sqrt[b^2 - 2*c*(a + c) + b*Sqrt[b^2 - 4*a*c]]) + (b*Cos[x])/c^2 - (Cos[x]*Sin[x])/(2*c)

________________________________________________________________________________________

Rubi [A]  time = 3.12929, antiderivative size = 323, normalized size of antiderivative = 1., number of steps used = 12, number of rules used = 8, integrand size = 19, \(\frac{\text{number of rules}}{\text{integrand size}}\) = 0.421, Rules used = {3256, 2638, 2635, 8, 3292, 2660, 618, 204} \[ -\frac{\sqrt{2} \left (-\frac{2 a^2 c^2-4 a b^2 c+b^4}{\sqrt{b^2-4 a c}}-2 a b c+b^3\right ) \tan ^{-1}\left (\frac{\tan \left (\frac{x}{2}\right ) \left (b-\sqrt{b^2-4 a c}\right )+2 c}{\sqrt{2} \sqrt{-b \sqrt{b^2-4 a c}-2 c (a+c)+b^2}}\right )}{c^3 \sqrt{-b \sqrt{b^2-4 a c}-2 c (a+c)+b^2}}-\frac{\sqrt{2} \left (\frac{2 a^2 c^2-4 a b^2 c+b^4}{\sqrt{b^2-4 a c}}-2 a b c+b^3\right ) \tan ^{-1}\left (\frac{\tan \left (\frac{x}{2}\right ) \left (\sqrt{b^2-4 a c}+b\right )+2 c}{\sqrt{2} \sqrt{b \sqrt{b^2-4 a c}-2 c (a+c)+b^2}}\right )}{c^3 \sqrt{b \sqrt{b^2-4 a c}-2 c (a+c)+b^2}}+\frac{x \left (b^2-a c\right )}{c^3}+\frac{b \cos (x)}{c^2}+\frac{x}{2 c}-\frac{\sin (x) \cos (x)}{2 c} \]

Antiderivative was successfully verified.

[In]

Int[Sin[x]^4/(a + b*Sin[x] + c*Sin[x]^2),x]

[Out]

x/(2*c) + ((b^2 - a*c)*x)/c^3 - (Sqrt[2]*(b^3 - 2*a*b*c - (b^4 - 4*a*b^2*c + 2*a^2*c^2)/Sqrt[b^2 - 4*a*c])*Arc
Tan[(2*c + (b - Sqrt[b^2 - 4*a*c])*Tan[x/2])/(Sqrt[2]*Sqrt[b^2 - 2*c*(a + c) - b*Sqrt[b^2 - 4*a*c]])])/(c^3*Sq
rt[b^2 - 2*c*(a + c) - b*Sqrt[b^2 - 4*a*c]]) - (Sqrt[2]*(b^3 - 2*a*b*c + (b^4 - 4*a*b^2*c + 2*a^2*c^2)/Sqrt[b^
2 - 4*a*c])*ArcTan[(2*c + (b + Sqrt[b^2 - 4*a*c])*Tan[x/2])/(Sqrt[2]*Sqrt[b^2 - 2*c*(a + c) + b*Sqrt[b^2 - 4*a
*c]])])/(c^3*Sqrt[b^2 - 2*c*(a + c) + b*Sqrt[b^2 - 4*a*c]]) + (b*Cos[x])/c^2 - (Cos[x]*Sin[x])/(2*c)

Rule 3256

Int[sin[(d_.) + (e_.)*(x_)]^(m_.)*((a_.) + (b_.)*sin[(d_.) + (e_.)*(x_)]^(n_.) + (c_.)*sin[(d_.) + (e_.)*(x_)]
^(n2_.))^(p_), x_Symbol] :> Int[ExpandTrig[sin[d + e*x]^m*(a + b*sin[d + e*x]^n + c*sin[d + e*x]^(2*n))^p, x],
x] /; FreeQ[{a, b, c, d, e}, x] && EqQ[n2, 2*n] && NeQ[b^2 - 4*a*c, 0] && IntegersQ[m, n, p]

Rule 2638

Int[sin[(c_.) + (d_.)*(x_)], x_Symbol] :> -Simp[Cos[c + d*x]/d, x] /; FreeQ[{c, d}, x]

Rule 2635

Int[((b_.)*sin[(c_.) + (d_.)*(x_)])^(n_), x_Symbol] :> -Simp[(b*Cos[c + d*x]*(b*Sin[c + d*x])^(n - 1))/(d*n),
x] + Dist[(b^2*(n - 1))/n, Int[(b*Sin[c + d*x])^(n - 2), x], x] /; FreeQ[{b, c, d}, x] && GtQ[n, 1] && Integer
Q[2*n]

Rule 8

Int[a_, x_Symbol] :> Simp[a*x, x] /; FreeQ[a, x]

Rule 3292

Int[((A_) + (B_.)*sin[(d_.) + (e_.)*(x_)])/((a_.) + (b_.)*sin[(d_.) + (e_.)*(x_)] + (c_.)*sin[(d_.) + (e_.)*(x
_)]^2), x_Symbol] :> Module[{q = Rt[b^2 - 4*a*c, 2]}, Dist[B + (b*B - 2*A*c)/q, Int[1/(b + q + 2*c*Sin[d + e*x
]), x], x] + Dist[B - (b*B - 2*A*c)/q, Int[1/(b - q + 2*c*Sin[d + e*x]), x], x]] /; FreeQ[{a, b, c, d, e, A, B
}, x] && NeQ[b^2 - 4*a*c, 0]

Rule 2660

Int[((a_) + (b_.)*sin[(c_.) + (d_.)*(x_)])^(-1), x_Symbol] :> With[{e = FreeFactors[Tan[(c + d*x)/2], x]}, Dis
t[(2*e)/d, Subst[Int[1/(a + 2*b*e*x + a*e^2*x^2), x], x, Tan[(c + d*x)/2]/e], x]] /; FreeQ[{a, b, c, d}, x] &&
NeQ[a^2 - b^2, 0]

Rule 618

Int[((a_.) + (b_.)*(x_) + (c_.)*(x_)^2)^(-1), x_Symbol] :> Dist[-2, Subst[Int[1/Simp[b^2 - 4*a*c - x^2, x], x]
, x, b + 2*c*x], x] /; FreeQ[{a, b, c}, x] && NeQ[b^2 - 4*a*c, 0]

Rule 204

Int[((a_) + (b_.)*(x_)^2)^(-1), x_Symbol] :> -Simp[ArcTan[(Rt[-b, 2]*x)/Rt[-a, 2]]/(Rt[-a, 2]*Rt[-b, 2]), x] /
; FreeQ[{a, b}, x] && PosQ[a/b] && (LtQ[a, 0] || LtQ[b, 0])

Rubi steps

\begin{align*} \int \frac{\sin ^4(x)}{a+b \sin (x)+c \sin ^2(x)} \, dx &=\int \left (\frac{b^2-a c}{c^3}-\frac{b \sin (x)}{c^2}+\frac{\sin ^2(x)}{c}+\frac{-a b^2 \left (1-\frac{a c}{b^2}\right )-b^3 \left (1-\frac{2 a c}{b^2}\right ) \sin (x)}{c^3 \left (a+b \sin (x)+c \sin ^2(x)\right )}\right ) \, dx\\ &=\frac{\left (b^2-a c\right ) x}{c^3}+\frac{\int \frac{-a b^2 \left (1-\frac{a c}{b^2}\right )-b^3 \left (1-\frac{2 a c}{b^2}\right ) \sin (x)}{a+b \sin (x)+c \sin ^2(x)} \, dx}{c^3}-\frac{b \int \sin (x) \, dx}{c^2}+\frac{\int \sin ^2(x) \, dx}{c}\\ &=\frac{\left (b^2-a c\right ) x}{c^3}+\frac{b \cos (x)}{c^2}-\frac{\cos (x) \sin (x)}{2 c}+\frac{\int 1 \, dx}{2 c}-\frac{\left (b^3-2 a b c-\frac{b^4-4 a b^2 c+2 a^2 c^2}{\sqrt{b^2-4 a c}}\right ) \int \frac{1}{b-\sqrt{b^2-4 a c}+2 c \sin (x)} \, dx}{c^3}-\frac{\left (b^3-2 a b c+\frac{b^4-4 a b^2 c+2 a^2 c^2}{\sqrt{b^2-4 a c}}\right ) \int \frac{1}{b+\sqrt{b^2-4 a c}+2 c \sin (x)} \, dx}{c^3}\\ &=\frac{x}{2 c}+\frac{\left (b^2-a c\right ) x}{c^3}+\frac{b \cos (x)}{c^2}-\frac{\cos (x) \sin (x)}{2 c}-\frac{\left (2 \left (b^3-2 a b c-\frac{b^4-4 a b^2 c+2 a^2 c^2}{\sqrt{b^2-4 a c}}\right )\right ) \operatorname{Subst}\left (\int \frac{1}{b-\sqrt{b^2-4 a c}+4 c x+\left (b-\sqrt{b^2-4 a c}\right ) x^2} \, dx,x,\tan \left (\frac{x}{2}\right )\right )}{c^3}-\frac{\left (2 \left (b^3-2 a b c+\frac{b^4-4 a b^2 c+2 a^2 c^2}{\sqrt{b^2-4 a c}}\right )\right ) \operatorname{Subst}\left (\int \frac{1}{b+\sqrt{b^2-4 a c}+4 c x+\left (b+\sqrt{b^2-4 a c}\right ) x^2} \, dx,x,\tan \left (\frac{x}{2}\right )\right )}{c^3}\\ &=\frac{x}{2 c}+\frac{\left (b^2-a c\right ) x}{c^3}+\frac{b \cos (x)}{c^2}-\frac{\cos (x) \sin (x)}{2 c}+\frac{\left (4 \left (b^3-2 a b c-\frac{b^4-4 a b^2 c+2 a^2 c^2}{\sqrt{b^2-4 a c}}\right )\right ) \operatorname{Subst}\left (\int \frac{1}{-8 \left (b^2-2 c (a+c)-b \sqrt{b^2-4 a c}\right )-x^2} \, dx,x,4 c+2 \left (b-\sqrt{b^2-4 a c}\right ) \tan \left (\frac{x}{2}\right )\right )}{c^3}+\frac{\left (4 \left (b^3-2 a b c+\frac{b^4-4 a b^2 c+2 a^2 c^2}{\sqrt{b^2-4 a c}}\right )\right ) \operatorname{Subst}\left (\int \frac{1}{4 \left (4 c^2-\left (b+\sqrt{b^2-4 a c}\right )^2\right )-x^2} \, dx,x,4 c+2 \left (b+\sqrt{b^2-4 a c}\right ) \tan \left (\frac{x}{2}\right )\right )}{c^3}\\ &=\frac{x}{2 c}+\frac{\left (b^2-a c\right ) x}{c^3}-\frac{\sqrt{2} \left (b^3-2 a b c-\frac{b^4-4 a b^2 c+2 a^2 c^2}{\sqrt{b^2-4 a c}}\right ) \tan ^{-1}\left (\frac{2 c+\left (b-\sqrt{b^2-4 a c}\right ) \tan \left (\frac{x}{2}\right )}{\sqrt{2} \sqrt{b^2-2 c (a+c)-b \sqrt{b^2-4 a c}}}\right )}{c^3 \sqrt{b^2-2 c (a+c)-b \sqrt{b^2-4 a c}}}-\frac{\sqrt{2} \left (b^3-2 a b c+\frac{b^4-4 a b^2 c+2 a^2 c^2}{\sqrt{b^2-4 a c}}\right ) \tan ^{-1}\left (\frac{2 c+\left (b+\sqrt{b^2-4 a c}\right ) \tan \left (\frac{x}{2}\right )}{\sqrt{2} \sqrt{b^2-2 c (a+c)+b \sqrt{b^2-4 a c}}}\right )}{c^3 \sqrt{b^2-2 c (a+c)+b \sqrt{b^2-4 a c}}}+\frac{b \cos (x)}{c^2}-\frac{\cos (x) \sin (x)}{2 c}\\ \end{align*}

Mathematica [C]  time = 1.3118, size = 410, normalized size = 1.27 \[ \frac{-\frac{4 \left (2 i a^2 c^2+b^3 \sqrt{4 a c-b^2}-4 i a b^2 c-2 a b c \sqrt{4 a c-b^2}+i b^4\right ) \tan ^{-1}\left (\frac{2 c+\tan \left (\frac{x}{2}\right ) \left (b-i \sqrt{4 a c-b^2}\right )}{\sqrt{2} \sqrt{-i b \sqrt{4 a c-b^2}-2 c (a+c)+b^2}}\right )}{\sqrt{2 a c-\frac{b^2}{2}} \sqrt{-i b \sqrt{4 a c-b^2}-2 c (a+c)+b^2}}-\frac{4 \left (-2 i a^2 c^2+b^3 \sqrt{4 a c-b^2}+4 i a b^2 c-2 a b c \sqrt{4 a c-b^2}-i b^4\right ) \tan ^{-1}\left (\frac{2 c+\tan \left (\frac{x}{2}\right ) \left (b+i \sqrt{4 a c-b^2}\right )}{\sqrt{2} \sqrt{i b \sqrt{4 a c-b^2}-2 c (a+c)+b^2}}\right )}{\sqrt{2 a c-\frac{b^2}{2}} \sqrt{i b \sqrt{4 a c-b^2}-2 c (a+c)+b^2}}+2 c x (c-2 a)+4 b^2 x+4 b c \cos (x)-c^2 \sin (2 x)}{4 c^3} \]

Antiderivative was successfully verified.

[In]

Integrate[Sin[x]^4/(a + b*Sin[x] + c*Sin[x]^2),x]

[Out]

(4*b^2*x + 2*c*(-2*a + c)*x - (4*(I*b^4 - (4*I)*a*b^2*c + (2*I)*a^2*c^2 + b^3*Sqrt[-b^2 + 4*a*c] - 2*a*b*c*Sqr
t[-b^2 + 4*a*c])*ArcTan[(2*c + (b - I*Sqrt[-b^2 + 4*a*c])*Tan[x/2])/(Sqrt[2]*Sqrt[b^2 - 2*c*(a + c) - I*b*Sqrt
[-b^2 + 4*a*c]])])/(Sqrt[-b^2/2 + 2*a*c]*Sqrt[b^2 - 2*c*(a + c) - I*b*Sqrt[-b^2 + 4*a*c]]) - (4*((-I)*b^4 + (4
*I)*a*b^2*c - (2*I)*a^2*c^2 + b^3*Sqrt[-b^2 + 4*a*c] - 2*a*b*c*Sqrt[-b^2 + 4*a*c])*ArcTan[(2*c + (b + I*Sqrt[-
b^2 + 4*a*c])*Tan[x/2])/(Sqrt[2]*Sqrt[b^2 - 2*c*(a + c) + I*b*Sqrt[-b^2 + 4*a*c]])])/(Sqrt[-b^2/2 + 2*a*c]*Sqr
t[b^2 - 2*c*(a + c) + I*b*Sqrt[-b^2 + 4*a*c]]) + 4*b*c*Cos[x] - c^2*Sin[2*x])/(4*c^3)

________________________________________________________________________________________

Maple [B]  time = 0.188, size = 1181, normalized size = 3.7 \begin{align*} \text{result too large to display} \end{align*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

int(sin(x)^4/(a+b*sin(x)+c*sin(x)^2),x)

[Out]

12/c^2*a^2/(8*a*c-2*b^2)/(4*c*a-2*b^2+2*b*(-4*a*c+b^2)^(1/2)+4*a^2)^(1/2)*arctan((-2*a*tan(1/2*x)+(-4*a*c+b^2)
^(1/2)-b)/(4*c*a-2*b^2+2*b*(-4*a*c+b^2)^(1/2)+4*a^2)^(1/2))*(-4*a*c+b^2)^(1/2)*b-4/c^3*a/(8*a*c-2*b^2)/(4*c*a-
2*b^2+2*b*(-4*a*c+b^2)^(1/2)+4*a^2)^(1/2)*arctan((-2*a*tan(1/2*x)+(-4*a*c+b^2)^(1/2)-b)/(4*c*a-2*b^2+2*b*(-4*a
*c+b^2)^(1/2)+4*a^2)^(1/2))*(-4*a*c+b^2)^(1/2)*b^3-16/c*a^3/(8*a*c-2*b^2)/(4*c*a-2*b^2+2*b*(-4*a*c+b^2)^(1/2)+
4*a^2)^(1/2)*arctan((-2*a*tan(1/2*x)+(-4*a*c+b^2)^(1/2)-b)/(4*c*a-2*b^2+2*b*(-4*a*c+b^2)^(1/2)+4*a^2)^(1/2))+2
0/c^2*a^2/(8*a*c-2*b^2)/(4*c*a-2*b^2+2*b*(-4*a*c+b^2)^(1/2)+4*a^2)^(1/2)*arctan((-2*a*tan(1/2*x)+(-4*a*c+b^2)^
(1/2)-b)/(4*c*a-2*b^2+2*b*(-4*a*c+b^2)^(1/2)+4*a^2)^(1/2))*b^2-4/c^3*a/(8*a*c-2*b^2)/(4*c*a-2*b^2+2*b*(-4*a*c+
b^2)^(1/2)+4*a^2)^(1/2)*arctan((-2*a*tan(1/2*x)+(-4*a*c+b^2)^(1/2)-b)/(4*c*a-2*b^2+2*b*(-4*a*c+b^2)^(1/2)+4*a^
2)^(1/2))*b^4+12/c^2*a^2/(8*a*c-2*b^2)/(4*c*a-2*b^2-2*b*(-4*a*c+b^2)^(1/2)+4*a^2)^(1/2)*arctan((2*a*tan(1/2*x)
+b+(-4*a*c+b^2)^(1/2))/(4*c*a-2*b^2-2*b*(-4*a*c+b^2)^(1/2)+4*a^2)^(1/2))*(-4*a*c+b^2)^(1/2)*b-4/c^3*a/(8*a*c-2
*b^2)/(4*c*a-2*b^2-2*b*(-4*a*c+b^2)^(1/2)+4*a^2)^(1/2)*arctan((2*a*tan(1/2*x)+b+(-4*a*c+b^2)^(1/2))/(4*c*a-2*b
^2-2*b*(-4*a*c+b^2)^(1/2)+4*a^2)^(1/2))*(-4*a*c+b^2)^(1/2)*b^3+16/c*a^3/(8*a*c-2*b^2)/(4*c*a-2*b^2-2*b*(-4*a*c
+b^2)^(1/2)+4*a^2)^(1/2)*arctan((2*a*tan(1/2*x)+b+(-4*a*c+b^2)^(1/2))/(4*c*a-2*b^2-2*b*(-4*a*c+b^2)^(1/2)+4*a^
2)^(1/2))-20/c^2*a^2/(8*a*c-2*b^2)/(4*c*a-2*b^2-2*b*(-4*a*c+b^2)^(1/2)+4*a^2)^(1/2)*arctan((2*a*tan(1/2*x)+b+(
-4*a*c+b^2)^(1/2))/(4*c*a-2*b^2-2*b*(-4*a*c+b^2)^(1/2)+4*a^2)^(1/2))*b^2+4/c^3*a/(8*a*c-2*b^2)/(4*c*a-2*b^2-2*
b*(-4*a*c+b^2)^(1/2)+4*a^2)^(1/2)*arctan((2*a*tan(1/2*x)+b+(-4*a*c+b^2)^(1/2))/(4*c*a-2*b^2-2*b*(-4*a*c+b^2)^(
1/2)+4*a^2)^(1/2))*b^4+1/c/(tan(1/2*x)^2+1)^2*tan(1/2*x)^3+2/c^2/(tan(1/2*x)^2+1)^2*tan(1/2*x)^2*b-1/c/(tan(1/
2*x)^2+1)^2*tan(1/2*x)+2/c^2/(tan(1/2*x)^2+1)^2*b-2/c^2*arctan(tan(1/2*x))*a+2/c^3*arctan(tan(1/2*x))*b^2+1/2*
x/c

________________________________________________________________________________________

Maxima [F(-1)]  time = 0., size = 0, normalized size = 0. \begin{align*} \text{Timed out} \end{align*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate(sin(x)^4/(a+b*sin(x)+c*sin(x)^2),x, algorithm="maxima")

[Out]

Timed out

________________________________________________________________________________________

Fricas [B]  time = 36.3275, size = 16731, normalized size = 51.8 \begin{align*} \text{result too large to display} \end{align*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate(sin(x)^4/(a+b*sin(x)+c*sin(x)^2),x, algorithm="fricas")

[Out]

-1/4*(sqrt(2)*c^3*sqrt((a^2*b^6 - b^8 - 2*a^4*c^4 - 2*(a^5 - 8*a^3*b^2)*c^3 + (9*a^4*b^2 - 20*a^2*b^4)*c^2 - 2
*(3*a^3*b^4 - 4*a*b^6)*c + (4*a*c^9 + (8*a^2 - b^2)*c^8 + 2*(2*a^3 - 3*a*b^2)*c^7 - (a^2*b^2 - b^4)*c^6)*sqrt(
-(a^4*b^10 - 2*a^2*b^12 + b^14 + 16*a^6*b^2*c^6 + 8*(3*a^7*b^2 - 10*a^5*b^4)*c^5 + (9*a^8*b^2 - 92*a^6*b^4 + 1
48*a^4*b^6)*c^4 - 4*(6*a^7*b^4 - 31*a^5*b^6 + 32*a^3*b^8)*c^3 + 2*(11*a^6*b^6 - 37*a^4*b^8 + 28*a^2*b^10)*c^2
- 4*(2*a^5*b^8 - 5*a^3*b^10 + 3*a*b^12)*c)/(4*a*c^17 + (16*a^2 - b^2)*c^16 + 12*(2*a^3 - a*b^2)*c^15 + 2*(8*a^
4 - 11*a^2*b^2 + b^4)*c^14 + 4*(a^5 - 3*a^3*b^2 + 2*a*b^4)*c^13 - (a^4*b^2 - 2*a^2*b^4 + b^6)*c^12)))/(4*a*c^9
+ (8*a^2 - b^2)*c^8 + 2*(2*a^3 - 3*a*b^2)*c^7 - (a^2*b^2 - b^4)*c^6))*log(16*a^7*b*c^4 + 4*(3*a^8*b - 10*a^6*
b^3)*c^3 - 8*(2*a^7*b^3 - 3*a^5*b^5)*c^2 + 2*(4*a^5*c^9 + (8*a^6 - a^4*b^2)*c^8 + 2*(2*a^7 - 3*a^5*b^2)*c^7 -
(a^6*b^2 - a^4*b^4)*c^6)*sqrt(-(a^4*b^10 - 2*a^2*b^12 + b^14 + 16*a^6*b^2*c^6 + 8*(3*a^7*b^2 - 10*a^5*b^4)*c^5
+ (9*a^8*b^2 - 92*a^6*b^4 + 148*a^4*b^6)*c^4 - 4*(6*a^7*b^4 - 31*a^5*b^6 + 32*a^3*b^8)*c^3 + 2*(11*a^6*b^6 -
37*a^4*b^8 + 28*a^2*b^10)*c^2 - 4*(2*a^5*b^8 - 5*a^3*b^10 + 3*a*b^12)*c)/(4*a*c^17 + (16*a^2 - b^2)*c^16 + 12*
(2*a^3 - a*b^2)*c^15 + 2*(8*a^4 - 11*a^2*b^2 + b^4)*c^14 + 4*(a^5 - 3*a^3*b^2 + 2*a*b^4)*c^13 - (a^4*b^2 - 2*a
^2*b^4 + b^6)*c^12))*sin(x) + 4*(a^6*b^5 - a^4*b^7)*c - sqrt(2)*((8*a^3*c^12 + 6*(4*a^4 - 3*a^2*b^2)*c^11 + 2*
(12*a^5 - 25*a^3*b^2 + 4*a*b^4)*c^10 + (8*a^6 - 38*a^4*b^2 + 35*a^2*b^4 - b^6)*c^9 - 2*(3*a^5*b^2 - 8*a^3*b^4
+ 5*a*b^6)*c^8 + (a^4*b^4 - 2*a^2*b^6 + b^8)*c^7)*sqrt(-(a^4*b^10 - 2*a^2*b^12 + b^14 + 16*a^6*b^2*c^6 + 8*(3*
a^7*b^2 - 10*a^5*b^4)*c^5 + (9*a^8*b^2 - 92*a^6*b^4 + 148*a^4*b^6)*c^4 - 4*(6*a^7*b^4 - 31*a^5*b^6 + 32*a^3*b^
8)*c^3 + 2*(11*a^6*b^6 - 37*a^4*b^8 + 28*a^2*b^10)*c^2 - 4*(2*a^5*b^8 - 5*a^3*b^10 + 3*a*b^12)*c)/(4*a*c^17 +
(16*a^2 - b^2)*c^16 + 12*(2*a^3 - a*b^2)*c^15 + 2*(8*a^4 - 11*a^2*b^2 + b^4)*c^14 + 4*(a^5 - 3*a^3*b^2 + 2*a*b
^4)*c^13 - (a^4*b^2 - 2*a^2*b^4 + b^6)*c^12))*cos(x) - (32*a^5*b^2*c^6 + 8*(5*a^6*b^2 - 13*a^4*b^4)*c^5 + 2*(6
*a^7*b^2 - 47*a^5*b^4 + 56*a^3*b^6)*c^4 - (19*a^6*b^4 - 69*a^4*b^6 + 54*a^2*b^8)*c^3 + 4*(2*a^5*b^6 - 5*a^3*b^
8 + 3*a*b^10)*c^2 - (a^4*b^8 - 2*a^2*b^10 + b^12)*c)*cos(x))*sqrt((a^2*b^6 - b^8 - 2*a^4*c^4 - 2*(a^5 - 8*a^3*
b^2)*c^3 + (9*a^4*b^2 - 20*a^2*b^4)*c^2 - 2*(3*a^3*b^4 - 4*a*b^6)*c + (4*a*c^9 + (8*a^2 - b^2)*c^8 + 2*(2*a^3
- 3*a*b^2)*c^7 - (a^2*b^2 - b^4)*c^6)*sqrt(-(a^4*b^10 - 2*a^2*b^12 + b^14 + 16*a^6*b^2*c^6 + 8*(3*a^7*b^2 - 10
*a^5*b^4)*c^5 + (9*a^8*b^2 - 92*a^6*b^4 + 148*a^4*b^6)*c^4 - 4*(6*a^7*b^4 - 31*a^5*b^6 + 32*a^3*b^8)*c^3 + 2*(
11*a^6*b^6 - 37*a^4*b^8 + 28*a^2*b^10)*c^2 - 4*(2*a^5*b^8 - 5*a^3*b^10 + 3*a*b^12)*c)/(4*a*c^17 + (16*a^2 - b^
2)*c^16 + 12*(2*a^3 - a*b^2)*c^15 + 2*(8*a^4 - 11*a^2*b^2 + b^4)*c^14 + 4*(a^5 - 3*a^3*b^2 + 2*a*b^4)*c^13 - (
a^4*b^2 - 2*a^2*b^4 + b^6)*c^12)))/(4*a*c^9 + (8*a^2 - b^2)*c^8 + 2*(2*a^3 - 3*a*b^2)*c^7 - (a^2*b^2 - b^4)*c^
6)) + 2*(a^6*b^6 - a^4*b^8 + 4*a^7*b^2*c^3 + (3*a^8*b^2 - 10*a^6*b^4)*c^2 - 2*(2*a^7*b^4 - 3*a^5*b^6)*c)*sin(x
)) - sqrt(2)*c^3*sqrt((a^2*b^6 - b^8 - 2*a^4*c^4 - 2*(a^5 - 8*a^3*b^2)*c^3 + (9*a^4*b^2 - 20*a^2*b^4)*c^2 - 2*
(3*a^3*b^4 - 4*a*b^6)*c - (4*a*c^9 + (8*a^2 - b^2)*c^8 + 2*(2*a^3 - 3*a*b^2)*c^7 - (a^2*b^2 - b^4)*c^6)*sqrt(-
(a^4*b^10 - 2*a^2*b^12 + b^14 + 16*a^6*b^2*c^6 + 8*(3*a^7*b^2 - 10*a^5*b^4)*c^5 + (9*a^8*b^2 - 92*a^6*b^4 + 14
8*a^4*b^6)*c^4 - 4*(6*a^7*b^4 - 31*a^5*b^6 + 32*a^3*b^8)*c^3 + 2*(11*a^6*b^6 - 37*a^4*b^8 + 28*a^2*b^10)*c^2 -
4*(2*a^5*b^8 - 5*a^3*b^10 + 3*a*b^12)*c)/(4*a*c^17 + (16*a^2 - b^2)*c^16 + 12*(2*a^3 - a*b^2)*c^15 + 2*(8*a^4
- 11*a^2*b^2 + b^4)*c^14 + 4*(a^5 - 3*a^3*b^2 + 2*a*b^4)*c^13 - (a^4*b^2 - 2*a^2*b^4 + b^6)*c^12)))/(4*a*c^9
+ (8*a^2 - b^2)*c^8 + 2*(2*a^3 - 3*a*b^2)*c^7 - (a^2*b^2 - b^4)*c^6))*log(16*a^7*b*c^4 + 4*(3*a^8*b - 10*a^6*b
^3)*c^3 - 8*(2*a^7*b^3 - 3*a^5*b^5)*c^2 - 2*(4*a^5*c^9 + (8*a^6 - a^4*b^2)*c^8 + 2*(2*a^7 - 3*a^5*b^2)*c^7 - (
a^6*b^2 - a^4*b^4)*c^6)*sqrt(-(a^4*b^10 - 2*a^2*b^12 + b^14 + 16*a^6*b^2*c^6 + 8*(3*a^7*b^2 - 10*a^5*b^4)*c^5
+ (9*a^8*b^2 - 92*a^6*b^4 + 148*a^4*b^6)*c^4 - 4*(6*a^7*b^4 - 31*a^5*b^6 + 32*a^3*b^8)*c^3 + 2*(11*a^6*b^6 - 3
7*a^4*b^8 + 28*a^2*b^10)*c^2 - 4*(2*a^5*b^8 - 5*a^3*b^10 + 3*a*b^12)*c)/(4*a*c^17 + (16*a^2 - b^2)*c^16 + 12*(
2*a^3 - a*b^2)*c^15 + 2*(8*a^4 - 11*a^2*b^2 + b^4)*c^14 + 4*(a^5 - 3*a^3*b^2 + 2*a*b^4)*c^13 - (a^4*b^2 - 2*a^
2*b^4 + b^6)*c^12))*sin(x) + 4*(a^6*b^5 - a^4*b^7)*c - sqrt(2)*((8*a^3*c^12 + 6*(4*a^4 - 3*a^2*b^2)*c^11 + 2*(
12*a^5 - 25*a^3*b^2 + 4*a*b^4)*c^10 + (8*a^6 - 38*a^4*b^2 + 35*a^2*b^4 - b^6)*c^9 - 2*(3*a^5*b^2 - 8*a^3*b^4 +
5*a*b^6)*c^8 + (a^4*b^4 - 2*a^2*b^6 + b^8)*c^7)*sqrt(-(a^4*b^10 - 2*a^2*b^12 + b^14 + 16*a^6*b^2*c^6 + 8*(3*a
^7*b^2 - 10*a^5*b^4)*c^5 + (9*a^8*b^2 - 92*a^6*b^4 + 148*a^4*b^6)*c^4 - 4*(6*a^7*b^4 - 31*a^5*b^6 + 32*a^3*b^8
)*c^3 + 2*(11*a^6*b^6 - 37*a^4*b^8 + 28*a^2*b^10)*c^2 - 4*(2*a^5*b^8 - 5*a^3*b^10 + 3*a*b^12)*c)/(4*a*c^17 + (
16*a^2 - b^2)*c^16 + 12*(2*a^3 - a*b^2)*c^15 + 2*(8*a^4 - 11*a^2*b^2 + b^4)*c^14 + 4*(a^5 - 3*a^3*b^2 + 2*a*b^
4)*c^13 - (a^4*b^2 - 2*a^2*b^4 + b^6)*c^12))*cos(x) + (32*a^5*b^2*c^6 + 8*(5*a^6*b^2 - 13*a^4*b^4)*c^5 + 2*(6*
a^7*b^2 - 47*a^5*b^4 + 56*a^3*b^6)*c^4 - (19*a^6*b^4 - 69*a^4*b^6 + 54*a^2*b^8)*c^3 + 4*(2*a^5*b^6 - 5*a^3*b^8
+ 3*a*b^10)*c^2 - (a^4*b^8 - 2*a^2*b^10 + b^12)*c)*cos(x))*sqrt((a^2*b^6 - b^8 - 2*a^4*c^4 - 2*(a^5 - 8*a^3*b
^2)*c^3 + (9*a^4*b^2 - 20*a^2*b^4)*c^2 - 2*(3*a^3*b^4 - 4*a*b^6)*c - (4*a*c^9 + (8*a^2 - b^2)*c^8 + 2*(2*a^3 -
3*a*b^2)*c^7 - (a^2*b^2 - b^4)*c^6)*sqrt(-(a^4*b^10 - 2*a^2*b^12 + b^14 + 16*a^6*b^2*c^6 + 8*(3*a^7*b^2 - 10*
a^5*b^4)*c^5 + (9*a^8*b^2 - 92*a^6*b^4 + 148*a^4*b^6)*c^4 - 4*(6*a^7*b^4 - 31*a^5*b^6 + 32*a^3*b^8)*c^3 + 2*(1
1*a^6*b^6 - 37*a^4*b^8 + 28*a^2*b^10)*c^2 - 4*(2*a^5*b^8 - 5*a^3*b^10 + 3*a*b^12)*c)/(4*a*c^17 + (16*a^2 - b^2
)*c^16 + 12*(2*a^3 - a*b^2)*c^15 + 2*(8*a^4 - 11*a^2*b^2 + b^4)*c^14 + 4*(a^5 - 3*a^3*b^2 + 2*a*b^4)*c^13 - (a
^4*b^2 - 2*a^2*b^4 + b^6)*c^12)))/(4*a*c^9 + (8*a^2 - b^2)*c^8 + 2*(2*a^3 - 3*a*b^2)*c^7 - (a^2*b^2 - b^4)*c^6
)) + 2*(a^6*b^6 - a^4*b^8 + 4*a^7*b^2*c^3 + (3*a^8*b^2 - 10*a^6*b^4)*c^2 - 2*(2*a^7*b^4 - 3*a^5*b^6)*c)*sin(x)
) + sqrt(2)*c^3*sqrt((a^2*b^6 - b^8 - 2*a^4*c^4 - 2*(a^5 - 8*a^3*b^2)*c^3 + (9*a^4*b^2 - 20*a^2*b^4)*c^2 - 2*(
3*a^3*b^4 - 4*a*b^6)*c - (4*a*c^9 + (8*a^2 - b^2)*c^8 + 2*(2*a^3 - 3*a*b^2)*c^7 - (a^2*b^2 - b^4)*c^6)*sqrt(-(
a^4*b^10 - 2*a^2*b^12 + b^14 + 16*a^6*b^2*c^6 + 8*(3*a^7*b^2 - 10*a^5*b^4)*c^5 + (9*a^8*b^2 - 92*a^6*b^4 + 148
*a^4*b^6)*c^4 - 4*(6*a^7*b^4 - 31*a^5*b^6 + 32*a^3*b^8)*c^3 + 2*(11*a^6*b^6 - 37*a^4*b^8 + 28*a^2*b^10)*c^2 -
4*(2*a^5*b^8 - 5*a^3*b^10 + 3*a*b^12)*c)/(4*a*c^17 + (16*a^2 - b^2)*c^16 + 12*(2*a^3 - a*b^2)*c^15 + 2*(8*a^4
- 11*a^2*b^2 + b^4)*c^14 + 4*(a^5 - 3*a^3*b^2 + 2*a*b^4)*c^13 - (a^4*b^2 - 2*a^2*b^4 + b^6)*c^12)))/(4*a*c^9 +
(8*a^2 - b^2)*c^8 + 2*(2*a^3 - 3*a*b^2)*c^7 - (a^2*b^2 - b^4)*c^6))*log(-16*a^7*b*c^4 - 4*(3*a^8*b - 10*a^6*b
^3)*c^3 + 8*(2*a^7*b^3 - 3*a^5*b^5)*c^2 + 2*(4*a^5*c^9 + (8*a^6 - a^4*b^2)*c^8 + 2*(2*a^7 - 3*a^5*b^2)*c^7 - (
a^6*b^2 - a^4*b^4)*c^6)*sqrt(-(a^4*b^10 - 2*a^2*b^12 + b^14 + 16*a^6*b^2*c^6 + 8*(3*a^7*b^2 - 10*a^5*b^4)*c^5
+ (9*a^8*b^2 - 92*a^6*b^4 + 148*a^4*b^6)*c^4 - 4*(6*a^7*b^4 - 31*a^5*b^6 + 32*a^3*b^8)*c^3 + 2*(11*a^6*b^6 - 3
7*a^4*b^8 + 28*a^2*b^10)*c^2 - 4*(2*a^5*b^8 - 5*a^3*b^10 + 3*a*b^12)*c)/(4*a*c^17 + (16*a^2 - b^2)*c^16 + 12*(
2*a^3 - a*b^2)*c^15 + 2*(8*a^4 - 11*a^2*b^2 + b^4)*c^14 + 4*(a^5 - 3*a^3*b^2 + 2*a*b^4)*c^13 - (a^4*b^2 - 2*a^
2*b^4 + b^6)*c^12))*sin(x) - 4*(a^6*b^5 - a^4*b^7)*c - sqrt(2)*((8*a^3*c^12 + 6*(4*a^4 - 3*a^2*b^2)*c^11 + 2*(
12*a^5 - 25*a^3*b^2 + 4*a*b^4)*c^10 + (8*a^6 - 38*a^4*b^2 + 35*a^2*b^4 - b^6)*c^9 - 2*(3*a^5*b^2 - 8*a^3*b^4 +
5*a*b^6)*c^8 + (a^4*b^4 - 2*a^2*b^6 + b^8)*c^7)*sqrt(-(a^4*b^10 - 2*a^2*b^12 + b^14 + 16*a^6*b^2*c^6 + 8*(3*a
^7*b^2 - 10*a^5*b^4)*c^5 + (9*a^8*b^2 - 92*a^6*b^4 + 148*a^4*b^6)*c^4 - 4*(6*a^7*b^4 - 31*a^5*b^6 + 32*a^3*b^8
)*c^3 + 2*(11*a^6*b^6 - 37*a^4*b^8 + 28*a^2*b^10)*c^2 - 4*(2*a^5*b^8 - 5*a^3*b^10 + 3*a*b^12)*c)/(4*a*c^17 + (
16*a^2 - b^2)*c^16 + 12*(2*a^3 - a*b^2)*c^15 + 2*(8*a^4 - 11*a^2*b^2 + b^4)*c^14 + 4*(a^5 - 3*a^3*b^2 + 2*a*b^
4)*c^13 - (a^4*b^2 - 2*a^2*b^4 + b^6)*c^12))*cos(x) + (32*a^5*b^2*c^6 + 8*(5*a^6*b^2 - 13*a^4*b^4)*c^5 + 2*(6*
a^7*b^2 - 47*a^5*b^4 + 56*a^3*b^6)*c^4 - (19*a^6*b^4 - 69*a^4*b^6 + 54*a^2*b^8)*c^3 + 4*(2*a^5*b^6 - 5*a^3*b^8
+ 3*a*b^10)*c^2 - (a^4*b^8 - 2*a^2*b^10 + b^12)*c)*cos(x))*sqrt((a^2*b^6 - b^8 - 2*a^4*c^4 - 2*(a^5 - 8*a^3*b
^2)*c^3 + (9*a^4*b^2 - 20*a^2*b^4)*c^2 - 2*(3*a^3*b^4 - 4*a*b^6)*c - (4*a*c^9 + (8*a^2 - b^2)*c^8 + 2*(2*a^3 -
3*a*b^2)*c^7 - (a^2*b^2 - b^4)*c^6)*sqrt(-(a^4*b^10 - 2*a^2*b^12 + b^14 + 16*a^6*b^2*c^6 + 8*(3*a^7*b^2 - 10*
a^5*b^4)*c^5 + (9*a^8*b^2 - 92*a^6*b^4 + 148*a^4*b^6)*c^4 - 4*(6*a^7*b^4 - 31*a^5*b^6 + 32*a^3*b^8)*c^3 + 2*(1
1*a^6*b^6 - 37*a^4*b^8 + 28*a^2*b^10)*c^2 - 4*(2*a^5*b^8 - 5*a^3*b^10 + 3*a*b^12)*c)/(4*a*c^17 + (16*a^2 - b^2
)*c^16 + 12*(2*a^3 - a*b^2)*c^15 + 2*(8*a^4 - 11*a^2*b^2 + b^4)*c^14 + 4*(a^5 - 3*a^3*b^2 + 2*a*b^4)*c^13 - (a
^4*b^2 - 2*a^2*b^4 + b^6)*c^12)))/(4*a*c^9 + (8*a^2 - b^2)*c^8 + 2*(2*a^3 - 3*a*b^2)*c^7 - (a^2*b^2 - b^4)*c^6
)) - 2*(a^6*b^6 - a^4*b^8 + 4*a^7*b^2*c^3 + (3*a^8*b^2 - 10*a^6*b^4)*c^2 - 2*(2*a^7*b^4 - 3*a^5*b^6)*c)*sin(x)
) - sqrt(2)*c^3*sqrt((a^2*b^6 - b^8 - 2*a^4*c^4 - 2*(a^5 - 8*a^3*b^2)*c^3 + (9*a^4*b^2 - 20*a^2*b^4)*c^2 - 2*(
3*a^3*b^4 - 4*a*b^6)*c + (4*a*c^9 + (8*a^2 - b^2)*c^8 + 2*(2*a^3 - 3*a*b^2)*c^7 - (a^2*b^2 - b^4)*c^6)*sqrt(-(
a^4*b^10 - 2*a^2*b^12 + b^14 + 16*a^6*b^2*c^6 + 8*(3*a^7*b^2 - 10*a^5*b^4)*c^5 + (9*a^8*b^2 - 92*a^6*b^4 + 148
*a^4*b^6)*c^4 - 4*(6*a^7*b^4 - 31*a^5*b^6 + 32*a^3*b^8)*c^3 + 2*(11*a^6*b^6 - 37*a^4*b^8 + 28*a^2*b^10)*c^2 -
4*(2*a^5*b^8 - 5*a^3*b^10 + 3*a*b^12)*c)/(4*a*c^17 + (16*a^2 - b^2)*c^16 + 12*(2*a^3 - a*b^2)*c^15 + 2*(8*a^4
- 11*a^2*b^2 + b^4)*c^14 + 4*(a^5 - 3*a^3*b^2 + 2*a*b^4)*c^13 - (a^4*b^2 - 2*a^2*b^4 + b^6)*c^12)))/(4*a*c^9 +
(8*a^2 - b^2)*c^8 + 2*(2*a^3 - 3*a*b^2)*c^7 - (a^2*b^2 - b^4)*c^6))*log(-16*a^7*b*c^4 - 4*(3*a^8*b - 10*a^6*b
^3)*c^3 + 8*(2*a^7*b^3 - 3*a^5*b^5)*c^2 - 2*(4*a^5*c^9 + (8*a^6 - a^4*b^2)*c^8 + 2*(2*a^7 - 3*a^5*b^2)*c^7 - (
a^6*b^2 - a^4*b^4)*c^6)*sqrt(-(a^4*b^10 - 2*a^2*b^12 + b^14 + 16*a^6*b^2*c^6 + 8*(3*a^7*b^2 - 10*a^5*b^4)*c^5
+ (9*a^8*b^2 - 92*a^6*b^4 + 148*a^4*b^6)*c^4 - 4*(6*a^7*b^4 - 31*a^5*b^6 + 32*a^3*b^8)*c^3 + 2*(11*a^6*b^6 - 3
7*a^4*b^8 + 28*a^2*b^10)*c^2 - 4*(2*a^5*b^8 - 5*a^3*b^10 + 3*a*b^12)*c)/(4*a*c^17 + (16*a^2 - b^2)*c^16 + 12*(
2*a^3 - a*b^2)*c^15 + 2*(8*a^4 - 11*a^2*b^2 + b^4)*c^14 + 4*(a^5 - 3*a^3*b^2 + 2*a*b^4)*c^13 - (a^4*b^2 - 2*a^
2*b^4 + b^6)*c^12))*sin(x) - 4*(a^6*b^5 - a^4*b^7)*c - sqrt(2)*((8*a^3*c^12 + 6*(4*a^4 - 3*a^2*b^2)*c^11 + 2*(
12*a^5 - 25*a^3*b^2 + 4*a*b^4)*c^10 + (8*a^6 - 38*a^4*b^2 + 35*a^2*b^4 - b^6)*c^9 - 2*(3*a^5*b^2 - 8*a^3*b^4 +
5*a*b^6)*c^8 + (a^4*b^4 - 2*a^2*b^6 + b^8)*c^7)*sqrt(-(a^4*b^10 - 2*a^2*b^12 + b^14 + 16*a^6*b^2*c^6 + 8*(3*a
^7*b^2 - 10*a^5*b^4)*c^5 + (9*a^8*b^2 - 92*a^6*b^4 + 148*a^4*b^6)*c^4 - 4*(6*a^7*b^4 - 31*a^5*b^6 + 32*a^3*b^8
)*c^3 + 2*(11*a^6*b^6 - 37*a^4*b^8 + 28*a^2*b^10)*c^2 - 4*(2*a^5*b^8 - 5*a^3*b^10 + 3*a*b^12)*c)/(4*a*c^17 + (
16*a^2 - b^2)*c^16 + 12*(2*a^3 - a*b^2)*c^15 + 2*(8*a^4 - 11*a^2*b^2 + b^4)*c^14 + 4*(a^5 - 3*a^3*b^2 + 2*a*b^
4)*c^13 - (a^4*b^2 - 2*a^2*b^4 + b^6)*c^12))*cos(x) - (32*a^5*b^2*c^6 + 8*(5*a^6*b^2 - 13*a^4*b^4)*c^5 + 2*(6*
a^7*b^2 - 47*a^5*b^4 + 56*a^3*b^6)*c^4 - (19*a^6*b^4 - 69*a^4*b^6 + 54*a^2*b^8)*c^3 + 4*(2*a^5*b^6 - 5*a^3*b^8
+ 3*a*b^10)*c^2 - (a^4*b^8 - 2*a^2*b^10 + b^12)*c)*cos(x))*sqrt((a^2*b^6 - b^8 - 2*a^4*c^4 - 2*(a^5 - 8*a^3*b
^2)*c^3 + (9*a^4*b^2 - 20*a^2*b^4)*c^2 - 2*(3*a^3*b^4 - 4*a*b^6)*c + (4*a*c^9 + (8*a^2 - b^2)*c^8 + 2*(2*a^3 -
3*a*b^2)*c^7 - (a^2*b^2 - b^4)*c^6)*sqrt(-(a^4*b^10 - 2*a^2*b^12 + b^14 + 16*a^6*b^2*c^6 + 8*(3*a^7*b^2 - 10*
a^5*b^4)*c^5 + (9*a^8*b^2 - 92*a^6*b^4 + 148*a^4*b^6)*c^4 - 4*(6*a^7*b^4 - 31*a^5*b^6 + 32*a^3*b^8)*c^3 + 2*(1
1*a^6*b^6 - 37*a^4*b^8 + 28*a^2*b^10)*c^2 - 4*(2*a^5*b^8 - 5*a^3*b^10 + 3*a*b^12)*c)/(4*a*c^17 + (16*a^2 - b^2
)*c^16 + 12*(2*a^3 - a*b^2)*c^15 + 2*(8*a^4 - 11*a^2*b^2 + b^4)*c^14 + 4*(a^5 - 3*a^3*b^2 + 2*a*b^4)*c^13 - (a
^4*b^2 - 2*a^2*b^4 + b^6)*c^12)))/(4*a*c^9 + (8*a^2 - b^2)*c^8 + 2*(2*a^3 - 3*a*b^2)*c^7 - (a^2*b^2 - b^4)*c^6
)) - 2*(a^6*b^6 - a^4*b^8 + 4*a^7*b^2*c^3 + (3*a^8*b^2 - 10*a^6*b^4)*c^2 - 2*(2*a^7*b^4 - 3*a^5*b^6)*c)*sin(x)
) + 2*c^2*cos(x)*sin(x) - 4*b*c*cos(x) - 2*(2*b^2 - 2*a*c + c^2)*x)/c^3

________________________________________________________________________________________

Sympy [F(-1)]  time = 0., size = 0, normalized size = 0. \begin{align*} \text{Timed out} \end{align*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate(sin(x)**4/(a+b*sin(x)+c*sin(x)**2),x)

[Out]

Timed out

________________________________________________________________________________________

Giac [F(-1)]  time = 0., size = 0, normalized size = 0. \begin{align*} \text{Timed out} \end{align*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate(sin(x)^4/(a+b*sin(x)+c*sin(x)^2),x, algorithm="giac")

[Out]

Timed out